Non-adiabatic current generation in a finite width semiconductor ring 
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We consider a model of a semiconductor quantum ring of finite width in a constant perpendicular 
magnetic field. We show how a current of the same order as the persistent current can be generated 
non-adiabatically by a short intensive pulse in the Tera-Hertz regime. 



Time-dependent or radiation induced phenomena in 
various mesoscopic electron systems have caught the at- 
tention of experimental^ and theoretical^ groups recently, 
just to mention two cases of a long Hst. The generation of 
a current by a slowly oscillating harmonic magnetic flux 
in a strictly one-dimensional (ID) ballistic quantum ring 
with a spinless Luttinger liquid has been investigated by 
Moskalets^i finding that the AC part depends both on 
the driving frequency and the temperature, but the DC 
part of the current is only very sensitive to the driving 
amplitude. More recently Moskalets and Biittiker have 
discovered a strong parity effect in a ID quantum wire 
excited by a harmonic delta-function potential)^ They 
find the current in a ring with an odd number of spinless 
electrons to be diamagnetic, exactly like in a pure ring 
without an excitation, but in a ring with an even number 
the current oscillates in sign with a large amplitude and 
with a small period compared to the unit fiux quantum 
$0 = h/{ec). 

Inspired by the orbital current modes found in ellipti- 
cal quantum dots by Serra et sl.^S^ and the large effects 
of a nonlinear excitation of quantum dotsp we present 
here an alternative method to generate DC currents in 
quantum rings nonadiabatically with a strong and short 
THz pulse. In our case the finite width of the quantum 
ring is essential and the current generation does rely on 
a change in the many-electron structure of the system 
caused by the THz pulse. 

The equilibrium persistent currents and the magneti- 
zation of quantum rings of finite width has been studied 
by Tan and Inksonyi and Vasile et al.j^ and their low- 
intensity far-infrared absorption has been measured by 
Lorke et al.;^ and calculated by several groups piSiiiii^ 

To calculate the time evolution of the system requires 
the knowledge of the initial state, its ground state, at a 
particular time t — to, which within a suitable mean field 
approximation is described by 

H\a) = {Ho + H, + Vc + Hint) \a) = e^la), (1) 

where Hq is the Hamiltonian for the 2DEG, that is sub- 
ject to the perpendicular magnetic field B = Bz, and 
confined by the parabolic potential Vpni{r) — m*WQr^/2. 
We modify the radial confinement with a central poten- 



Hcr = {l/2)g* fiBBaz, where g* is an effective 5-factor, 
fiB the Bohr magneton, and az = ±1 for the two spin 
projections. 

The H amiltonian Hq defines a natural length scale, a — 
l^/uj^, in terms of the magnetic length, I = yJhc/{eB), 
the cyclotron frequency lOc — eB/(rn*c), and the charac- 
teristic frequency, SI — ^10^ + AbJ^ . The eigenstates |a), 
with a = (M, Ur), of the Fock Darwin^ Hamiltonian Hq 
are used as a mathematical basis in which the eigenstates 
of H, I a), are expanded in, and in which (QJ is trans- 
formed into a matrix eigenvalue equation. No restriction 
is put on the angular symmetry of the system. A suit- 
able truncation of the basis that can be varied to 
improve the accuracy, leads to a numerically tractable 
eigenvalue problem. 

kt t — to the quantum ring is radiated by a short 
THz pulse making the Hamiltonian of the system time- 
dependent, H{t) = H + W(t), with 

W{t) = Ftrl^pl cos (iVp0) exp (-sr^ - Vt) 

sin (wii) sin (ajt)0(7r — wi^), (2) 

where Q is the Heaviside step function. The parameters, 
s, lu, and Np shall be explained below. Since the 
radiation is of arbitrary strength we resume to the task 
of solving the equation of motion for the density operator 



[H + W{t),p{t)]. 



(3) 



The structure of this equation is inconvenient for numer- 
ical evaluation so we resort instead to the time-evolution 
operator T, defined by p{t) — T{t)poT^{t), which has 
the simpler equation of motion 



ihT{t) = H{t)T{t) 
-ihf+{t) = T+{t)H(t). 



(4) 



tial hill, Vc{r) 
potential V^onf 
ring like. The spin-dependent part of the Hamiltonian is 



— Vbexp(— 7r ), creating the confining 

— Vpar + Vc in order to make the system 



We discretize time and use the Crank-Nicholson algo- 
rithm for the time-integration with the initial condition, 
T(0) = 1. This is performed in the truncated Fock- 
Darwin basis {|a)} for the corresponding matrix version 
of the Eqs. lHJ. 

As a matter of convenience, we use the time-dependent 
orbital magnetization M.o(t) to quantify the AC and DC 
currents induced by W{t). The magnetization operator 
Mo = (l/2c)r X j together with the definition for the 
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current, j = — er = {—ie/h)[H{t),r], give in terms of 
the density matrix in the Fock-Darwin basis {\a)} the 
magnetization 

{{S\x\f3)mH{m{l\y\r) - (/3|y|7)(7|i?W|T)) 
-(5|y|/3)((/3|i/W|7)(7k|r) - (/3|x|7)(7|i?W|T))X5) 

In this communication, for reasons of clarity and space, 
we only report on results for non- interacting electrons. 
The inclusion of the Coulomb interaction in a LSDA (Lo- 
cal Spin-Density Approximation) manner only changes 
the results here quantitatively, and shall be reported else- 
where. It should though be stated here, since it is also 
relevant to the noninteracting model, that the accuracy 
and the fine-tuning of the LSDA model have been tested 
for the ground state with comparison to the exact results 
calculated for a circular parabolically confined quantum 
dotji^ and in the time-evolution by convincing us that the 
model can reproduce the Kohn results for the absorption 
of a Far-infrared radiation 

In order to model a ring-shaped 2DEG with a number 
of electrons, N — 1 — 13, in a GaAs sample with vanishing 
electron density in its center we select the parameters, 
hujo = 3.37 meV, Vq — 30 meV, and a^7 = 1.0 in the 
center part of the confinement potential Vc- We assume 
g = 0.44, and m* = 0.067m. As this small effective mass 
enhances the orbital magnetization of the 2DEG with a 
factor of 14.9 with respect to the spin magnetization we 
shall ignore the latter here. Moreover, we are here merely 
using the orbital magnetization of the system as a mea- 
sure of the persistent currents in it. For the radiation 
pulse ^ of just over 3 ps, we select F = 2 ps~^. The en- 
velope frequency corresponds to huji = 0.658 meV, and 
the base frequency is huj = 2.63 meV. A conventional 
dipole radiation is described by \Np\ = 1 in Eq. 
but we shall also consider here higher multipoles with 
|iVp| = 2,3, or the monopole, Np — 0, which could only 
be produced by specially patterned microwave antennae 
or gate fingers. For TVp = we use sa^ = 1.0, but oth- 
erwise s = 0. We use a basis set with M'^^^ = 12, and 
n™'^'^ e [4,8]. 

We integrate the equations for the time-evolution op- 
erator 1^1 for a time interval of 40 ps using an increment 
in the range of 0.002 - 0.01 ps and evaluate Aio{t) at each 
step. The Fourier transform of the time series in the in- 
terval 4 - 40 ps, when the external excitation has been 
turned off, yields information about the average magne- 
tization in this interval, (Aio), and the oscillations of 
the time-dependent part. The external pulse ^ always 
produces a radial force on the 2DEG starting radial os- 
cillations of its density in the confining potential Konf- 
For all the values of Np considered here this oscillation is 
characterized by one or few narrow resonances close to 7 
meV, even though the driving frequency was hoj = 2.63 
meV. In an interacting system this frequency is modified 
by the Coulomb interaction. In Fig. we show the equi- 



librium density for 12 electrons in a ring together with 
the induced density at a later time for Np = 3. On top 




FIG. 1: The electron density (left), and the induced density 
(right) at some later point in time. B — 0.6 T, |A'p| — 3, 
Vta'^"' = 1.0 meV, iV = 12, and T = 1 K. 

of the radial oscillation there is a clear angular oscilla- 
tion visible in Fig. The average magnetization (Mo) 
after the excitation pulse has been turned off is differ- 
ent from the magnetization of the corresponding ground 
state Mo{0) as shown in Fig.|2l Moreover, pulse strength 
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FIG. 2: The equilibrium magnetization A4o{0), {+), and 
the averaged time-dependent change in the magnetization 
AA4o = {Mo{t) — Mo{0)) after the radiation pulse has van- 
ished for \Np\ = 1 (□), liVpl = 2 (*), and \Np\ = 3 (x). 
For \Np\ = 2 and 3 Ka'^pl = 1.0 meV, and for \Np\ = 1 
Ftal^pl = 2.0 meV. For |iVp| = 1, and \Np\ = 2 the curves 
have been magnified by a factor of 10. TV = 12, T = 1 K, and 
other parameters are in text. 

here is large enough to increase the average size of the 
system, just as has been observed in highly excited quan- 
tum dots 4. We are in the nonlinear regime of excitation 
and have nonadiabatically pumped energy into the sys- 
tem. As a result the 2DEG is now in an excited state 
with a different persistent current. This conclusion is 
supported by comparing the density matrix before and 
after the excitation. 
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We can neither generate a nonzero current this way 
at vanishing magnetic field, B = 0, nor with a pure ra- 
dial excitation, iVp = at a finite magnetic field. The 
reason for the \Np\ = 3 excitation to be more effective 
here than the \Np\ — 2 excitation is partly due to their 
different radial strength, and partly due to the stronger 
angular modulation for the |A^p| = 3 case. The |A'^p = 1 
excitation is different in the sense that it also causes a 
large oscillation of the center of mass of the system, a 
dipole oscillation absent in the other cases. The change 
Vt —Vt is ineffective since it can always be canceled, 
for symmetry reasons, by a fixed-angle rotation of the 
2DEG, the numerical model also reproduces this result. 

The current generated depends strongly on the number 
of electrons in the ring as can be seen by comparing the 
results for one electron in Fig.Elto the ring with 12 elec- 
trons, Fig.[2| A certain parity effect can be observed in 
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FIG. 3: The equilibrium magnetization M-o{Q), (+), and the 
averaged time-dependent change in the magnetization AAlo 
after the radiation pulse has vanished for |A'p| = 3 (x). 
V^iftl'^f I = 1.0 meV, A/' = 1, and r = 1 K. 



the sign of A^o(O), or the direction of the induced persis- 
tent current, see Fig.^ The magnetization in the ground 
state A^o(O) assumes a maximum value for an even num- 
ber of noninteracting electrons with spin and is much 
smaller for an odd number. This behavior can be un- 
derstood in terms of the quasi periodicity of the energy 
spectrum for a low magnetic field, see Fig. |3^b). Each 
energy level is almost spin degenerate and, for example, 
the change from 6 to 8 electrons changes the system from 
diamagnetic to paramagnetic as can be guessed from the 
slope of the highest occupied level. In case of 7 electrons 
the 2DEG pertains both qualities resulting in a small net 
magnetization. The induced magnetization (A^o), on the 
other hand is in antiphase with these oscillations. This 
can only be understood in terms of a changed character 
of the excited 2DEG. If the excitation is not too large 
then the system gains through it more of the character 
of the lowest unoccupied state, that fulfills the selection 
rule |AM| = \Np\. Thus, a still stronger excitation might 
be able to reverse this trend, as can be confirmed in Fig. 
El where an increase in the excitation strength changes 
the 2DEG from diamagnetic to paramagnetic. 
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FIG. 4: The equilibrium magnetization 7Vfo(0), (+), and the 
averaged time-dependent change in the magnetization IS.M0 
after the radiation pulse has vanished for |A'j,| = 3 (x) as 
functions of the number of electrons A*'. Vto}^"^^ = 1.0 meV, 
B = 0.4 T, and r = 1 K. 
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FIG. 5: The noninteracting single-electron energy spectrum 
for a quantum dot, Vo = meV (a), and a quantum ring with 
Vo = 30 meV (b). Arrows connect some low energy states 
with the same values of angular quantum number M . T — 1 
K, see other parameters in text. 



By inspecting the diagonal elements of the ground 
state density matrix for the case of TV = 12, and B = 0.4 
T it can be seen that there is a strong overlap with the 
5 basis states with = 0, and M ~ —2, —1, +1, +2, +3 
and a partial overlap with (M, rt^) — (0,0) and (0,1) 
refiecting the center potential hill in the system, that 
pushes the low lying M = states to higher energy as 
can be seen in Fig. With the center potential hill 
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FIG. 6: The time-dependent change in the magnetization 
AA4 after the radiation pulse has vanished for |A'p| — 3 (+) 
as a function of the pulse strength Vt. T = 1 K. 

present the states can not be assigned the radial quan- 
tum number Ur anymore. After the system has been 
excited with a Vjo'^ — 1 meV pulse the overlapping with 
the states (3,0), (0,0), (0,2), (0,1), and (1,1) continu- 
ously changes the most, reflecting a strong radial oscil- 
lation with a smaller angular mode. The excited state 
does thus include M-states in a different ratio causing a 
changed persistent current pattern. The occupation of 
still higher M-states needs still stronger excitation. 

Although the change in the average persistent current 
is a general feature of the strong excitation described 
here, clearly the antiphase found in Fig. 31 is not a gen- 
eral feature of the system, but rather the relation between 
the Mo and (Mo) depends strongly on the relevant selec- 
tion rule, I AM I = |A'^p|, and the energy spectrum, which 
depends strongly on the magnetic flux and is not simply 
periodic in the flux for a ring of flnite width, see Fig. 
Etb). 

We have demonstrated here that the persistent cur- 
rent, or its manifestation in the time-dependent magne- 
tization, of a quantum ring of a finite width can be con- 
trolled by a strong external excitation pulse by infiicting 
a nonadiabatic change upon the many-electron state of 



the system. This current generation or change in the 
persistent current can thus not be described by a linear 
response theory. Due to the fast response of the system, 
at the level of excitation applied here, the change in its 
state occurs almost entirely during the time the excita- 
tion pulse is applied to the 2DEG. The average magne- 
tization (Mo), describing the new persistent current in 
the system is evaluated after these switch-on effects have 
died down. Eventually the energy pumped into the sys- 
tem by the excitation pulse will dissipate. We make no 
attempts to describe the dissipation process and various 
possible channels thereof in the present model, but an- 
ticipate the lifetime of the excited state - the collective 
mode - to be long enough to allow for several cycles of 
oscillations before at least the energy is carried away by 
electromagnetic radiation. 

We have tested the results communicated here within 
a Density-Functional Theory approach using the LSDA 
and found similar qualitative results for the limited values 
of the system parameters tested. It should be emphasized 
that the neglect of the Coulomb interaction here is done 
consistently in the ground and the excited state by inte- 
grating directly the equation of motion, something that 
is not easy to perform within a linear response theory 
without risiking a violation of some conservation laws. 
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